arXiv:astro-ph/0512019vl 1 Dec 2005 


Astronomy Sz Astrophysics manuscript no. paper @ ESO 2008 

February 5, 2008 


A population synthesis study on the faint X-ray sources in the 

Galactic center region 

X.-W. Liu and X.-D. Li 

Department of Astronomy, Nanjing University, Nanjing 210093, China 
e-mail: liuxw@nju.edu.cn, lixd@nju.edu.cn 

Received March 29, 2005; accepted November 10, 2005 

ABSTRACT 

Context. Recent Chandra observations of the Galactic center region (GCR) have uncovered a population of 
sources. A few theoretical works have been made to investigate the nature of these sources. 

Aims. We examine the contributions and luminosity functions of various kinds of candidate objects which 
by previous authors or by ourselves. 

Methods. We conduct a population synthesis calculation based on Hurley et al.’s rapid binary evolution code. Several candidate 
models, i.e. wind-accreting neutron stars, intermediate polars, low mass X-ray binaries, young pulsars and massive stars with 
strong winds, are incorporated into our calculation. We also take the geometric effect of the accretion disk into account for 
Roche lobe overflow X-ray binaries. 

Results. Our results show that neutron star low-mass X-ray binaries contribute significantly to the observed sources. 
We also point out that wind-accreting neutron stars contribute negligibly to these sources due to propeller effect, and 
the intermediate polars play a minor role in accounting for the faint X-ray sources in both Wang et al. and Muno et al. 
survey. It should be mentioned that the majority of the sources in the survey field of Wang et al. are still beyond our expectation. 


faint discrete X-ray 
are proposed either 
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1. Introduction 

Instrumental development always leads people to study 
sources looking fainter. Recent Chandra observations of 
the Galactic center region (GCR) have uncovered a pop¬ 
ulation of faint discrete X-ray sources. Wang, Gotthelf, & 
Lang have reported hundreds of X-ray sources with 

luminosities Lx ~ 10^^ — 10^® erg s“^ in the GCR which 
contains ^ 1% of the total Galactic population (Pfahl 
et al. hereafter PRP02). Many of these sources are 

emitting in a hard energy (> 3 keV) band since the X-ray 
absorption, varying across the field, is particularly high 
in the 1 — 3 keV band. Muno et al. Il2()()3ll presented a 
catalog of 2357 point sources detected during their 590 ks 
Chandra observations of the 17' x 17' field around Sgr A*. 
It was estimated that there are ~ 2000 individual point¬ 
like sources with luminosities Lx ^ 10^° — 10^^ erg s“^ 
lying within a cylinder of radius 20 pc and depth 440 pc 
centered on the Galactic center at a distance of 8.5 kpc 
(Muno et al. I2004all . The spectra of these sources can be 
fitted by an absorbed power law with photon index P. 
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More than 1000 sources have relatively hard spectra with 

P < 1. 

The nature of the faint X-ray sources in the GGR has 
been investigated by several authors. Pfahl et al. II2()()2II at¬ 
tributed a significant fraction of the X-ray sources revealed 
by Wang et al. H2002II to wind-accreting neutron stars 
(WNSs). By calculating the flux distribution of WNSs, 
they concluded that the detected number ranges from ten 
to several hundred. The similar idea for pre-low-mass X- 
ray binaries was suggested by Willems & Kolb cm see 
however, PoDOv [^nM|l . Muno et al. H2004all proposed that 
magnetized cataclysmic variables (CVs), referred to as 
polars and intermediate polars (IPs), could account for 
^ 1000 — 2000 X-ray sources in their survey. The popula¬ 
tion synthesis work by Belczynski & Taam H2004L hereafter 
BT04) suggested that transient neutron star low-mass X- 
ray binaries (NS LMXBs) may contribute primarily to the 
X-ray source population in Muno et al. survey. 

In this paper, we employed the evolutionary popula¬ 
tion synthesis (EPS) method to calculate the expected 
numbers and luminosity distributions of various types of 
candidates that shine X-rays in Wang et al. and Muno 
et al. surveys. We examined the spin evolution of WNSs 
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proposed by PRP02, including the ejector, propeller and 
accretor stages (Davies & Pringle FlDHlI hereafter DP81; 
Lipunov Ikhsanov EoniiEnni- Besides traditional 

X-ray binaries and CVs, we also considered the contri¬ 
bution from rotation-powered pulsars and massive stars 
with strong winds. We describe the population synthesis 
method and the input physics for various types of X-ray 
sources in our model in §2. The calculated results are pre¬ 
sented in §3. Our discussions and conclusions are in §4. 


2. Model description 


We have used the EPS code developed by Hurley et al. 

to calculate the expected numbers in the 
Galactic disk for various types of single and binary X- 
ray source populations. This code incorporates evolution 
of single stars with binary-star interactions, such as mass 
transfer, mass accretion, common-envelope (CE) evolu¬ 
tion, collisions, supernova kicks, tidal friction and angular 
momentum loss mechanics. Most of our adopted parame¬ 
ters are the same as those described in Hurley et al. GlSi. 
We assume all stars are born in binary systems. The initial 
mass function (IMF) of Kroupa, Tout, & Gilmore lllhbdil is 
taken for the primary’s mass (Mi) distribution (discussed 
below in §2.3). For the secondary stars (M 2 ), we assume a 
uniform distribution of the mass ratio M 2 /Mi between 0 
and 1. A uniform distribution of In a is also taken for the 
binary separation a. The star formation rate parameter is 
S = 7.6085yr“^, corresponding to a rate of 0.02yr“^ 
for core-collapse supernovae (SNe) in our Galaxy, assume 
all the stars with masses > 8Mq die through SNe. During 
the SN explosions, a kick velocity Uk is imparted on the 
newborn compact stars with the Maxwellian distribution 





( 1 ) 


where a = 265 kms“^ (Hobbs et al. 120051) or 

190 kms“^ (Hansen & Phinney I1997II . The GE pa¬ 
rameter acE is set to be 1, if not otherwise mentioned. 
The velocity vw of winds from massive stars is taken to 
be Vw = ^/2P{GM/R), where the value of [3 depends 
on the spectral type of the stars. Hurley et al. cnna 
suggested it to be in the range of 0.125-7.0. 


2.1. Physical Models for X-ray Sources 

2.1.1. WNSs 

Similar as in PRP02, in our model a WNS is a neutron 
star (NS) accreting wind material from an intermediate- 
or high-mass (> 3Mq) un-evolved companion. We as¬ 
sume that the wind is steady and spherically symmet¬ 
ric. The standard Bondi-Hoyle accretion formula is 

employed to calculate the mass transfer rate to the NS. 
Be/NS binaries are not considered in the present model, 
due to their high eccentricities and the unclear structure 
of the Be star’s equatorial winds. PRP02 simply assumed 
that a WNS always accretes from the wind of its binary 


companion, ignoring its previous evolution as an ejector 
and propeller. Obviously the expected number of WNSs 
would be less than in PRP02, depending on the duration 
of the latter stages. The following equations give the esti¬ 
mates of the timescales for the ejector (a), supersonic (c) 
and subsonic (d) propeller phases for a natal NS in a wind 
environment (DP81; lkhsa,nov l20()2[l . 

Ta ~ 4.8 X 10®/r§'/M)/^/^/45E7“^ yr, (2) 

r, ~ 10 V§'i^MiV^"/ 45E7-' yr, (3) 

Td ~ 10^^;jfm/45Riooyr, (4) 

where ^ = lO^^/xai Gcm^ is the magnetic moment, P = 
lOOPioo s the spin period, Mns = wMq the mass, I = 
10 "*®/45 gcm^ the moment of inertia of the NS, respectively. 
V = lO'^Vr cms“^ is the relative velocity between the NS 
and the wind. The mass capture rate of the NS M = 
lO^’^Mir gs“^ is given by M = irR^pooV, where i?a = 
2 GMns/ 14^ is the accretion radius and poo is the density 
of the surrounding gas. When the spin period increases to 
the brake period Pbr = the NS 

switches from subsonic propeller to accretor. 

In the ejector stage, the transferred wind material is 
outside the light cylinder. The X-ray emission is generated 
from the rotational power of the NS. This will be discussed 
in §2.1.4. When the pressure of the ejected particles from 
the neutron star can no longer balance the ram pressure 
of the wind gas, the gas will penetrate into the light cylin¬ 
der and interact with the star’s magnetosphere, giving rise 
to X-ray emission. This propeller stage can be further di¬ 
vided into supersonic and subsonic stages (DP81). The 
corresponding X-ray luminosities are evaluated according 
to the work by Ikhsanov iiwmi . 

Finally when the magnetic pressure is less than the 
ram pressure of the falling material and the material at 
the magnetospheric boundary is cool enough so that the 
interchange instabilities can be triggered, the NS enters 
the accretor stage, and all the material captured by the NS 
is assumed to reach its surface. We then use Lx = pMc^ 
with ?7 ~ 0.1 to calculate the X-ray luminosity. 

The parameters of the natal NSs are set as follows. The 
initial spin periods P and magnetic fields B are chosen so 
that log P and log B are distributed normally with a mean 
of —2.3 and 12.5 respectively, and a standard deviation of 
0.3. 

2.1.2. Magnetic CVs 

Muno et al. Il2004all suggested that the majority of the 
GG sources are IPs. To examine this idea, we include GVs 
in our calculations, and set 5% of them are IPs (Kube et 

ai. i^uua . 

2.1.3. LMXBs 

This type of systems contain NSs or black holes (BHs) 
accreting from low-mass companions that overfill their 
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Roche lobes. Most of LMXBs are transients, with the 
majority of their time spent in quiescence with Lx < 
1034ej-gs“i. The criterion suggested by van Paradijs 
lllbbtill is used to determine whether the X-ray source is 
persistent or a soft X-ray transient (SXT). BT04 took 
a semi-empirical approach for the X-ray luminosities of 
LMXB transients in quiescence. We have instead adopted 
the theoretical model of Menou et al. for quies¬ 

cent SXTs, in which the inner part of the accretion flow is 
advection-dominated (ADAF, see Narayan, Mahadevan & 
Quataert 119981 for a review), surrounded by an outer thin 
disk. For BH SXTs, the quiescent luminosities are 

Lx,q = 5bh/ (1 + 5bh) X 10“® X ryMc^, (5) 

where gen = Madaf/M^^cc 0.5, Madaf is the rate at 
which mass is accreted via the ADAF, and Mace the rate 
at which mass is accumulated in the outer thin disk. The 
accumulated mass will be accreted in the outbursts. 

Menou et al. II1999II included the propeller effect to 
estimate the quiescent X-ray luminosities of NS transients 

Lx,q = gNs/(l + 9ns) ^ face X rjMc^, (6) 


2.1.5. Massive stars with strong winds 

Massive stars with strong winds, like Wolf-Rayet stars, 
are long known to be bright X-ray sources (Pallavicini 
et al. irrmni . The X-ray emission is generally thought to 
be generated in the lower layers of the winds, where the 
plasma is heated in shocks arising from small-scale wind 
structures that grow out from line-driven wind instabil¬ 
ities (Feldmeier et al. 119971 Dessart & Owocki 1200311 . 
Chelbowski & Garmany H1991I hereafter CG91) found 
that the X-ray luminosities of both Wolf-Rayet stars and 
O type stars can be fitted by the following relation, 

log Lx (erg s“^) = 34.32-|-0.34 log (Mwi^oo), (8) 

where Mw is the wind mass loss rate in units of M 0 yr“^ 
and Voc is the terminal velocity of the mass outflows in 
units of kms“^. 

Massive binaries can add considerable X-rays due to 
collision of the winds from the two stars under specific con¬ 
ditions. First, this X-ray enhancement is significant only 
when the binary separation a lies in the following range 
(CG91), 


where gxs — 0.2 and /acc is the fraction of the total mass 
through ADAF that reaches the NS surface. 

Both normal star and white dwarf (WD) donors are 
considered in our simulations. Apart from core collapse 
SNe, BT04 suggested the formation of NSs via accretion- 
induced collapse (AIG) of massive WDs. The EPS code in 
our simulation also includes this possibility: steady trans¬ 
fer of He-rich or C-rich material onto an ONe WD can lead 
to an AIC leaving an NS remnant, when the WD masses 
exceed the Chandrasekhar mass limit. 

2.1.4. Rotation-powered pulsars 

In the ejector stage, the spin-down of NSs is governed 
by the canonical rotation-powered pulsar mechanism. The 
rotation power is spent to generate the magnetic dipole 
waves and accelerate particles. Young pulsars are X-ray 
emitters. There appears to be a strong correlation between 
the rate of rotational energy loss E and their X-ray lumi¬ 
nosities Lx (Seward & Wang 119881 Becker & Trumper 
fTMTl Saito I1998II . The comprehensive investigation by 
Possenti et al. II2()()2II suggested that the X-ray luminosi¬ 
ties in 2 — 10 keV band depend on E with the following 
relation 

Lx = 10-i5-3(L;/erg erg s-h (7) 

Since E decreases with time, rotation-powered pulsars can 
shine in X-rays with luminosities ~ 10^^ — lO^^erg s“^ 
at ages of ~ 10^ — 10'^ yr. The X-ray flux contains both 
pulsed non-thermal X-ray emission from the NS magne¬ 
tosphere and unpulsed component from the synchrotron 
nebulae powered by the relativistic particles and magnetic 
helds ejected by the NS. Contribution from the NS cooling 
is most likely in the energy band < 2 keV and could be 
highly attenuated by the Galactic interstellar absorption. 


l.l(Ri -f i?2) < a < 5.5(Ri -f i?2), 


(9) 


where R is the stellar radius, and the subscripts 1 and 
2 denote the primary and secondary stars, respectively. 
Second, Prilutskii & Usov showed that the X-ray 

emission due to wind collision in binary systems is inef¬ 
ficient unless M 2 V 2 /{MiVi + M 2 V 2 ) > 0.4. If the above 
requirements are satisfied, we use the empirical relation 
of Portegies Zwart et al. Il2992afl to calculate the X-ray 
luminosities for colliding winds. 


Lx = 1.3 X 10^^( — 
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2.2. X-ray Band and Luminosities 

Wang et al. II2()()2II survey revealed ^ 1,000 discrete 
sources, most of which were detected in the relative hard 
energy range of 2-10 keV. Less than 20 of these sources are 
previously known (bright X-ray binaries). They suggested 
that half of these hard X-ray sources could be background 
active galactic nuclei (AGN). Considering the steep den¬ 
sity prohle in the Galactic disk, at least 100 sources with 
luminosities 10^^ - lO^^erg s“^ (Pfahl et al. 120021) are in the 
GCR. Among the 2357 point sources presented by Muno 
et al. imt . 281 sources were detected below 1.5 keV and 
are believed to be mainly in the foreground of the GC, and 
^100 sources are suggested to be background AGN. The 
surface density of the remaining sources, and their absorp¬ 
tion column toward the GC, in agreement with those in 
infrared surveys toward Sgr A*, demonstrate that these 
X-ray sources trace the general stellar population at the 
GCR (Muno et al. 1200311 . A tiny fraction of the remaining 
sources would be contaminated by the forground objects 
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because of the steep density profile of the Galactic disk. 
In fact, the ratio of the forground and the GCR sources 
(8.5 kpc X l)/(400pc X 1000) w 0.02 (see §2.4 for the 
GCR stellar density prescription), is negligible. By using 
the typical spectral model obtained from those relatively 
bright sources in this sample, Muno et al. I|200,*-{|l inferred 
the 2.0-8.0 keV luminosity range of 10^° ~ lO^^erg s“^ 
for the GC sources. The luminosities estimated for some 
of the models described in §2.1 are bolometric. They have 
to be converted into the 2-10 keV luminosities to be com¬ 
pared with observations. 

WNSs often exhibit non-thermal spectra that can be 
described with a T ~ 0 power law below 10 keV (e.g. 
Campana et al. 120011 ), we thus adopt their 2-10 keV lu¬ 
minosities to be half the calculated bolometric luminosi¬ 
ties. Magnetized CVs, especially IPs, usually exhibit hard 
thermal spectra {kT >10 keV) (e.g. Ezuka & lsbida lTT)!10|l . 
We also take their bolometric correction factor to be 0.5. 
The theoretical model of Menou et al. IIIOOOII of quies¬ 
cent X-ray transients is established for comparison with 
observational luminosities in 0.5-10 keV. So we multiply 
them by a factor of 0.5 to estimate the luminosities in 2- 
10 keV. The X-ray luminosities determined by Eq.(7) for 
rotation-powered pulsars are appropriately in 2-10 keV , 
therefore no correction is needed. GG91 obtained the em¬ 
pirical relation (Eq.8) for massive stars with strong winds 
from the data of Einstein observatory (in 0.5-3.5 keV). 
For comparison, Portegies Zwart et al. H2002a|l calculated 
the luminosities in the 0.5-3.5 keV band (Eq.lO) based on 
the archival Chandra X-ray observations. Taking account 
of their soft spectra {kT ~ 1 — 3 keV, Portegies Zwart et 
al. I2002a'|) . we choose the correction factor to be 0.1 for 
their 2-10 keV luminosities. 

The X-ray emission is generally assumed to be 
isotropic. However, for disk accreting sources, the geo¬ 
metric effect may affect the apparent luminosity distri¬ 
butions (Zhang, 12003) • Quantitatively, the intrinsic lumi¬ 
nosity is reduced by a factor of cos0(l-|-2cos0)/3, where 0 
is the inclination angle of the accretion disk; the factor of 
cos 9 is due to the area-projection effect, and the factor of 
(1-1-2 cos 9 )/3 is due to the limb-darkening effect (Netzer, 
[MB- The probability of seeing an accretion disk at an in¬ 
clination angle 9 is proportional to sin 9 if accretion disks 
are assumed to be oriented randomly in the sky. Therefore, 
the convolution between/(cc) = (1—x^)^/^[l-|-2(l— 

(where x = sin0 is uniformly distributed between 0 and 1) 
and a given intrinsic luminosity distribution produces the 
observed distribution (see e.g. Figures 2, 5 and 6). Thus, 
we take this inclination effect into account for Roche lobe 
overflow (RLOF) X-ray binaries (i.e. IP and LMXBs). 

2.3. Star Formation and IMF in GCR 

As noted by Muno et al. Il2004a|) . it is still a matter of 
debate as to whether the star formation is continuous or 
episodic, and whether it occurs only in localized regions 
or is relatively uniform throughout the GC. By modeling 


the evolution of the population of luminous infrared stars, 
Figer et al. mm suggested that the stat formation is 
probably continuous over the last ~ 10 Gyr, as we did in 
this paper. The Kroupa et al. II1993I) IMF adopted in this 
work has a slope of —2.7 when M > 1.0 Mq, derived from 
the stellar distribution towards both Galactic poles as well 
as the distribution of stars within 5.2 pc from the Sun. 
However, the molecular clouds in the GCR have higher 
densities and temperatures than those in the disk due to 
the tidal forces. This environment might favor the for¬ 
mation of massive stars (Morris 111221). The Hubble Space 
Telescope observation of Arches Cluster (which is within 
Wang et al. field) suggested an fairly flat IMF with a slope 
of ^ —1.7-—1.9 (Figer et al. 119991 Figer [21121 see also 
Stolte et al. mm . Portegies Zwart et al. (I2nn2hl) . how¬ 
ever, showed that the observed characteristics (unusually 
flat mass function and overabundance of massive stars) of 
the Arches cluster are consistent with a perfectly normal 
IMF. The observed anomalies are then caused by a combi¬ 
nation of observational selection effects and the dynamical 
evolution of the cluster. Considering lack of the knowledge 
of the star formation history in GCR, our assumptions of 
star formation and IMF in the GCR are not incompatible 
with recent observational and theoretical works. 


2.4. Expected Numbers 

The expected numbers for each class of X-ray sources 
in the GCR are related with the the ratio of the star 
population in the GCR and in the total disk. The size 
of Wang et al. survey held is 2° x 0.8°. Supposing 
the space density of stars in the Galactic disk to 
be n{R,z) ocexp(—i?/i?o)exp(—|z|/zo) and integrating 
n{R, z) over 1.6 deg^ through the GCR, Pfahl et al. II2()()2I) 
found that the held contains <1% of the total Galactic 
disk population. Alternatively, we can make an order of 
magnitude estimate on this ratio. Wang et al. survey 
held encompasses a physical area of 300 pc by 120 pc 
at a distance of 8.5 kpc. Assume that all of the point 
sources are within the nuclear bulge (which is about 300 
pc across in the radial direction; Mezger et al., 119961) and 
the average stellar density is 100 Mqpc“^ (which is 1000 
times that in the local neighborhood, Binney & Merriheld 
119981) . the Wang et al. held contains a stellar mass of 
300 X 120x300 X 1OO=1O®M0. Thus, this held contains 1% 
of the total Galactic stellar mass, which is ~ IO^^Mq. 
With EPS calculations, we try to search various types of 
X-ray sources (described in Sec.2) that can be detected in 
Wang et al. cnna and Muno et al. mm surveys. After 
we get the total number A^cai of a specihc type of sources 
in the Galactic disk, then ^ A^cai/lOO and ^ A^cai/lOO/d 
sources of this population are expected to be found in 
Wang et al. and Muno et al. surveys, respectively, since 
the latter held contains stars 4 times less than the former 
(BT04). Note that sources found in these two surveys have 
diherent luminosity ranges. 
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3. Results 

We adopt a variety of models (see Table 1), each with 
different assumptions for the parameters that govern the 
evolutions in the calculations. Tables 2 and 3 summarize 
the calculated numbers of various classes of X-ray sources 
contributed to Wang et al. and Muno et al. surveys, re¬ 
spectively. The expected numbers after correcting for the 
inclination effect are also listed. 

The luminosity distributions of WNSs are shown in 
Figure 1. The solid and dashed lines represent the calcu¬ 
lated results with and without the propeller effect con¬ 
sidered, respectively. The latter distribution implies that 
there are about 12 and 16 WNSs (for Models A and D 
respectively) in Wang et al. field. PRP02 predicted ~ 250 
WNSs in the solid angle of Wang et al field for their stan¬ 
dard kick model, and ~ 12 — 100 of them could be detected 
in the survey. This is compatible with our result (we get 
25 WNSs in Model F, in which the SN kick model is same 
as the standard model of PRP02). The figure also reveals 
that with decreasing wind velocities the WNS luminosities 
without considering the propeller effect increase, as ex¬ 
pected with the Bondi-Hoyle mechanism. However, when 
we take account of the propeller effect, the expected num¬ 
ber of WNSs that can he detected in Wang et al. survey is 
less than one, indicating that mass accretion may be highly 
inefficient during most part of the WNS stage. 

Figure 2 shows the luminosity distributions of vari¬ 
ous classes of RLOF X-ray binaries expected in Muno 
et al. field. The thin solid lines, dashed lines and thick 
solid lines are for IPs, NS LMXBs and BH LMXBs re¬ 
spectively. In this figure we can see that NS LMXBs in 
quiescence, generally brighter than BH LMXBs, account 
for a considerable fraction of the sources in Wang et al. 
survey (30 out of ~ 100 for model F, see Table 2); NS 
LMXBs with WD donor stars contribure a significant frac¬ 
tion of the sources in Muno et al. field. Most of the NS-WD 
LMXBs are ultracompact, with orbital periods less than 
2 hr. Their numbers (see Table 3) are roughly in line with 
that (291) in BT04. Our simulations predict hundreds of 
IPs in the GCR, fairly less than Muno et al. Il2004a|l ’s 
estimate (~ 1000). 

The CE efficiency parameter acE plays an important 
role in the evolution of close binaries. To see the depen¬ 
dence of the results on acE^ we can compare the results of 
Model A {acE = 1-0) with those of Model B {acE = 0.5) 
in Table 3. The parameter acE has two contrary effects on 
the evolution of close binaries. On one hand, high value of 
acE results in relatively wide binaries after the CE phase, 
so that the following RLOF mass transfer may be hard to 
take place. On the other hand, for the ultracompact X-ray 
binaries, a high value of a can prevent coalescence during 
the previous CE phase, significantly increasing the forma¬ 
tion rate of these systems. The overall influence is that 
larger a leads to more LMXBs, especially LMXBs with a 
WD donor, because these binaries have experienced CE 
evolution at least twice. 


Both IP and LMXB numbers increase significantly 
when the tidal spin-orbit coupling is considered in Models 
E and G. When the primary star ascends to the giant 
branch, its expansion causes a reduction of the orbital 
separation through the spin-orbit interaction, while wind 
mass loss from the evolved donor star causes the orbit to 
widen. The competition between the wind loss and the 
tidal interactions will determines the subsequent orbital 
evolution (Tauris & Savoniie l2()()()l Tauris l^ITITTIl . As a typ¬ 
ical example of CV formation in Model E (TIDE ON) let 
us consider binary stars of masses 2.9 and 0.4 Mq in an or¬ 
bit with initial separation of 484 Rq. About 0.5 Gyr later, 
the more massive star evolves to the AGB stage when 
the separation is 488 R©. In the following 3Myr evolu¬ 
tion , the tidal spin-orbit interaction dominates the orbital 
evolution and the binary separation decreases by 30% to 
355 Rq when the primary star begins to fill its Roche lobe. 
A CE forms and a CO white dwarf separated from the 
0.4 Mq main sequence (MS) star by 2Rq emerges. When 
the separation has fallen to 1.1 Rq (due to tidal interac¬ 
tion) 70Myr later, the secondary star fills its Roche lobe 
and the stage of CV evolution begins. In contrast, the 
same system in Model A (TIDE OFF) has an ever in¬ 
creasing orbital separation (from 484 to 500 Rq , owing to 
the wind mass loss effect) before the primary star fills its 
Roche lobe. The CE evolution leads to a CO WD-MS star 
binary with a separation of 3 Rq . The secondary star can¬ 
not fill its Roche lobe within 12 Gyr because the timescale 
for either its nuclear expansion or the orbital shrinkage 
due to gravitational radiation is too long. 

In Figure 3 we show the X-ray luminosity distributions 
of rotation-powered pulsars. There are totally ~28 and ~9 
pulsars predicted in Muno et al. and Wang et al. surveys 
respectively. The numbers of rotation-powered pulsars rely 
only on the core-collapse SN rate and their X-ray lifetime, 
both are independent of the parameters in Table 1. 

Figure 4 shows the X-ray luminosity distributions of 
massive stars with strong winds. Only a small fraction of 
them have colliding winds and thus are more luminous. 
Their luminosities increases with decreasing wind veloci¬ 
ties. The reason is that the X-ray luminosity from colliding 
winds is proportional to the product of the emissivity per 
unit volume, n^A, both the density n and emission rate A 
are anticorrelated with the wind velocity vw- We obtain 
^ 70 massive stars in Muno et al. survey. The number 
hardly changes with the model parameters. 

Figures 5 and 6 show the X-ray luminosity distri¬ 
butions of the total X-ray sources for different models. 
All models (especially models E and G with TIDE ON) 
seems to produce X-ray sources with reasonable numbers 
in Muno et al. field, but less than half the objects in Wang 
et al. field are accounted for. In Muno et al. the lu¬ 

minosity distribution peaks around 2 x lO^^erg s“^, nearly 
one order of magnitude higher than in our work. The rea¬ 
sons for this discrepancy can be addressed briefly as fol¬ 
lows. First the calculated mass transfer rates are long¬ 
term, averaged ones. They may not be directly compared 
with observed, instantaneous X-ray luminosities. Second, 
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observational selection effects may favor relatively more 
luminous X-ray sources to be identified. Third, some po¬ 
tentially important effects (such as the effect of X-ray ir¬ 
radiation on the secondary) have not been included in 
our binary calculations, which could considerably alter the 
evolution of these systems and increase the mass accretion 
rates (see also Podsiadlowski, Rappaport & Pfahl 120(1^ . 

We also calculate the numbers of bright X-ray sources 
in GCR. The outburst correction factor is chosen to be 
0.1 and 1 for the short- and long-period LMXB transient 
systems, respectively (see BT04 for details). We assume 
the bolometric correction factor is 0.1 and the duty cy¬ 
cle is 1%. For Model A, 23 X-ray sources brighter than 
10^®erg s“^ are found in Muno et al. field, including 16 
transient LMXBs in outburst and 7 persistent LMXBs. 
Only a few bright transients are detected in Muno et al. 
survey. For the Wang et al. field, we predict 92 bright 
point sources, whereas <20 bright sources were detected 
in Wang et al. II2002[I . Like BT04, our simulation seems to 
lead to an overproduction of the bright systems. However, 
for model B, the expected numbers of bright sources in 
Muno et al. and Wang et al. surveys are reduced to 9 and 
36, respectively. 

We calculate the point source contribution to the X- 
ray emission that has previously been ascribed to diffuse 
emission (Koyama et al. 119961 Sidoli & Mereghetti I1999II 
in the GCR. For various models, the X-ray point sources 
(< 10^®erg s“^) produce a mean surface brightness of 1.3— 
2.6 X 10“^^erg cm“^ s“^ arcmin“^. This is 3%-9% of that 
of the diffuse emission from the inner regions of the Galaxy 
derived by Sidoli & Mereghetti 1I1999II and by Koyama et 
al. II1996II . Our result is compatible with the observational 
result (10%) of Muno et al. H2()()3|l . 

4. Discussion and conclusions 

By use of the EPS method, we have investigated the na¬ 
ture of the faint X-ray sources in the GGR. Our sim¬ 
ulated objects include WNSs, magnetic GVs, LMXBs, 
rotation-powered pulsars and massive stars with strong 
winds. The main results are summarized in Tables 2 and 
3, with ^10-30 and ~1000-2000 sources expected to be 
detected in Wang et al. and Muno et al. surveys, respec¬ 
tively. For Wang et al. survey, we find that a consid¬ 
erable fraction of the discrete sources may be rotation- 
powered pulsars and NS LMXB transients in quiescence, 
while WNSs proposed by PRP02 have negligible contri¬ 
bution due to the propeller effect (Note that we have not 
considered Be/X-ray binaries in which the NS evolution 
may be quite different from investigated here.). Recent 
near infrared imaging of the X-ray sources in the GGR 
shows that the colour distribution of the identified candi¬ 
date counterparts of some of the X-ray sources is redder 
than expected for WNS systems, but consistent with later- 
type stars (Bandyopadhyay et al. I2()()5|l . For Muno et al. 
field, IPs present a minor contribution, and the majority 
of the X-ray sources seem to be NS LMXB transients with 
WD donors. The latter result is consistent with BT04. 


Our calculations also suggest that some of the point 
sources detected by Muno et al. survey may be massive 
stars with strong winds. Radio emission from these winds 
should be detectable at centimeter wavelengths (Panagia 
& Fell! 1197^1 Wright & Barlow fl975ll . We propose that 
a systematic radio observing campaign be undertaken to 
search for the stellar counterparts in the GGR surveyed by 
Muno et al. (some efforts have been made, e.g. by Lang 
et al. most recently). Moreover, radio observations 

could reveal the existence of jets, which are common in 
X-ray binaries. 

Sakano et al. recently reported the discov¬ 

ery in the GGR of two unusual X-ray transients XMM 
J174457-2850.3 and XMM J174544-2913.0 with flux 
variations in excess of a factor of 100 during roughly a 
year and peak X-ray luminosities of ~ 5 x 10^^ erg s“^ and 
with peculiar spectral features. Since no known classes of 
sources can well explain all their characteristics, these au¬ 
thors argued that these two sources may represent a new 
type of sources with different properties from those we 
have known. This also implies that what is happening in 
the GCR is much more complicated than we have already 
learnt. 
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Log (^'■9 s') 


Fig. 1. The X-ray luminosity distributions for WNSs with 
and without considering the propeller effect. Three Models 
C, A and D with decreasing wind velocity are shown from 
top to bottom. 
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Fig. 2. The X-ray luminosity distributions for Roche- 
lobe overflow systems including IPs (thin solid lines), NS 
LMXBs (dashed lines) and BH LMXBs (thick solid lines). 

The figure is plotted only for Model A since the main 
features do not change much as the parameters change. 

Upper panel: intrinsic distributions. Lower panel: distri¬ 
butions after correcting for accretion disk inclination ef¬ 
fect. 











10 


Liu & Li: Faint X-ray sources in GCR 



30 31 32 33 34 35 36 

Log (erg s ’) 


Fig. 3. The X-ray luminosity distributions for rotation- 
powered pulsars. This hgure is plot for Model A. 
Distributions for other models have little changes. 
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Fig. 4. The X-ray luminosity distributions for massive 
stars with strong winds (dashed lines). Distributions for 
wind-colliding binaries are highlighted with thick lines. 
Three Models C, A and D with decreasing wind velocity 
are shown from left to right. 
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Fig. 5. The total number of X-ray sources as a function 
of the X-ray luminosity. Three models, A (solid lines), B 
(dashed lines) and F (dotted lines), in which no tidal effect 
considered are shown. Upper panel: intrinsic distributions. 
Lower panel: distributions after correcting for the accre¬ 
tion disk inclination effect. 
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Fig. 6. Same as Fig.5. Two models, E (solid lines) and 
G (dashed lines), in which tidal effect was considered are 
shown here. Upper panel: intrinsic distributions. Lower 
panel: distributions after correcting for the accretion disk 
inclination effect. 
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Table 1. Model parameters. 


Model 

CecE 

/3 

a 

tides 

A 

1.0 

1.0 

265 

OFF 

B 

0.5 

1.0 

265 

OFF 

C 

1.0 

4.0 

265 

OFF 

D 

1.0 

0.25 

265 

OFF 

E 

1.0 

1.0 

265 

ON 

F 

1.0 

1.0 

190 

OFF 

G 

1.0 

1.0 

190 

ON 


Table 2. The expected numbers of X-ray sources in Wang 
et al. survey with Lx,2-iokeV = 10^^ — 10^®ergs“^. In 
parentheses we list the expected numbers of X-ray binaries 
with WD donors. The expected numbers after correcting 
for the inclination effect are listed in the last column. 


Model 

WNSs 

NS LMXBs 

BH LMXBs 

IPs 

Pulsars 

TOTAL 

TOTAL-inc 

A 

< 1 

3(0) 

<1(0) 

1 

9 

13 

13 

B 

< 1 

<1(0) 

<1(0) 

1 

9 

9 

9 

C 

< 1 

4(0) 

<1(0) 

1 

9 

14 

14 

D 

< 1 

<1(0) 

<1(0) 

2 

9 

11 

11 

E 

< 1 

2(0) 

1(0) 

1 

9 

13 

12 

F 

< 1 

23(0) 

<1(0) 

1 

9 

33 

30 

G 

< 1 

5(0) 

1(0) 

1 

9 

16 

14 


Table 3. The expected numbers of X-ray sources in Muno 
et al. survey with. Lx,2-iokeV = 10^° — 10^^ erg s“^. In 
parentheses we list the expected numbers of X-ray binaries 
with WD donors. The expected numbers after correcting 
for the inclination effect are listed in the last column. 


Model 

IPs 

NS LMXBs 

BH LMXBs 

Pulsars 

Massive Stars 

TOTAL 

TOTAL-inc 

A 

165 

648(646) 

132(0) 

28 

72 

1045 

912 

B 

169 

195(194) 

7(0) 

28 

71 

470 

405 

C 

163 

642(641) 

130(0) 

28 

72 

1035 

905 

D 

176 

935(934) 

158(0) 

28 

72 

1369 

1202 

E 

706 

1142(1141) 

270(0) 

28 

72 

2218 

1963 

F 

165 

694(693) 

128(0) 

28 

72 

1087 

961 

G 

706 

1242(1240) 

497(0) 

28 

71 

2544 

2248 





















